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Abstract Molecular phylogenetic study indicated that Luehdorfia chinensis huashanensis is ge- 
netically distinct from L. c. chinensis. Phylogenetic trees based on the COI and ND5+COI genes 
indicated that L. c. huashanensis is phylogenetically even more closely related to L. longicaudata 
than to L. c. chinensis. The present molecular phylogenetic study contributes a piece of evidence 
testifying to the significance of L. c. huashanensis and gives a basis for future conservation work. 
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Introduction 


Luehdorfia chinensis is widely distributed in central China, including Jiangsu, Zhejiang, 
Anwhei, Jiangxi, Henan, Hunan, Hubei, Shaanxi, and Sichuan provinces, and divided into 
two subspecies, L. c. huashanensis Lee and L. c. chinensis Leech (Watanabe 1996). L. c. 
huashanensis is distributed in the higher mountainous areas of the Hubei, Shaanxi and 
Sichuan provinces of China (Fig. 1) (Lee, 1982; Yuan et al., 1998) and L. c. chinensis is dis- 
tributed in the lower lands of Jiangsu, Zhejiang, Anwhei, Jiangxi, Henan, Hunan and Hubei 
provinces. Lee (1982) reported that morphological differences in these two subspecies 
were few in the larval and pupal stages, though L. c. huashanensis had distinct wing pat- 
terns, that is, narrower black and wider yellow bands in the adult (Watanabe, 1996). 
Furthermore, the two species feed on different plants; L. c. huashanensis on deciduous 
Asarum sieboldii, and L. c. chinensis on aiophyllous Asarum forbesii. The geographical 
distribution of L. c. huashanensis partially overlaps that of a closely related species, 
Luehdorfia longicaudata Lee. It is known, however, that these two species are well segre- 
gated by the fact that L. longicaudata feeds on Saruma henryi, and is distributed in lower 
altitude areas than L. c. huashanensis (Nose, 1990a; Watanabe, 1996). 


The speciation of the genus Luehdorfia has been studied from various viewpoints 
(Takahashi, 1973; Hiura, 1978; Ishizuka, 1980; Shinkawa, 1991). L. chinensis has also 
been characterized at the levels of male external genitalia (Shinkawa, 1991; Fujioka, 1992; 
Aoyama, 1994; Hong et al., 1999), adult antennae using SEM (Hong et al., 1997), cross- 
hybridization (Yashima et al., 1999) and DNA (Shinkawa, 1999; Makita et al., 2000). 
However, molecular studies on genetic differentiation among the subspecies of the genus 
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Fig. 1. Geographical distribution of Luehdorfia chinensis. 


Luehdorfia including L. c. huashanensis have not yet been reported. In the present study, 
phylogenetic analyses using mitochondrial DNA (mtDNA) are performed in order to char- 
acterize L. c. huashanensis. 


Materials and methods 


Samples 


Species of genus Luehdorfia collected for the present study are listed in Table 1 along with 
their collecting places. Sericinus montela and Papilio xuthus were used as out-group 
species for comparison. The nucleotide sequence data registered in the DDBJ data bank 
were also added for comparison. To prevent degradation of DNA in the specimens, samples 
were frozen under -20°C in zipped plastic bags immediately after squeezing the thorax. As 
for L. c. huashanensis, dried samples were employed. 


Methods 


MtDNA was extracted from the thoracic muscle of adults by using IsoQuick DNA extrac- 
tion kit (ORCA Research, Inc.). A389-bp region of the mitochondrial NADH dehydroge- 
nase subunit 5 (ND5) gene and a 457-bp of the cytochrome oxidase subunit I (COI) gene 
were amplified by PCR and the products obtained were treated with exonuclease I. Cycle 
sequencing reaction was carried out for the purified PCR products using the BigDye 
Terminator Cycle Sequence FS Ready Reaction Kit (Applied Biosystems) with the ABI 377 
DNA sequencer (Applied Biosystems). Primers V1 and C2 (Yagi et al., 1999) were used 
for PCR and direct sequencing of the NDS region. Similarly, the primers Ron and Nancy 
(Simon et al., 1994) were used for the COI region. Phylogenetic analysis of all samples 
were carried out using PAUP* 4.0b10 (Swoffoed, 1998). The neighbor-joining (NJ) method 
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Table 1. Luehdorfia and out-group species used in the present molecular phylogenetic study. 











DDBJ No. 











Species Code locality of collection NDS CO! References 
Luehdorfia japonica LI1 Tottori, Kokufu, Japan AB179879 AB179870 
LJ2 Yamaguchi, Kumage, Japan AB179880 AB179871 

Kanazawa, Japan — AF170867 
Sakyo-ku, Kyoto, Japan AB013142 — 

L. c. chinensis LCC Zhejiang, Hangzho, China AB016826 AB179872 

L. c. huashanesis LCH1 Shaanxi, Mt Qin Ling in Xi’an, China AB179881 AB179873 
Shaanxi, Mt Qin Ling in Xi’an, China AB179882 - 

L. longicudata LL Shaanxi, Mt Qin Ling in Xi’an, China AB016828 AB179874 

L. p. puziloi LPP Primorsky. Khabarovsk, Russia AB016829 AB179878 

L. p. coreana LPC Kyongsangnamdo. Nan Hae, Korea AB179883 AB179875 
Hwa San, Kyonggido, Korea — AB016830 

L. p. inexpecta LPI1 Nagano, Kamiina, Japan AB179884 AB179876 

LPI2 Aomori Hiraga, Japan AB179885 -— 

L. p. yessoensis LPY Hokkaido, Hamamasu, Japan AB179886 AB179877 

Allancastria cerisyi Outgroup AB095662 AF170869 1,4 

Bhutanitis mansfieldi Outgroup AB026727 — 3 

Papilio xuthus Outgroup AB013149 AF043999 2,6 

Parnassius apollo Outgroup AB095635 AY346254 4,5 

Sericinus montela Outgroup AB013144 AF170868 1,6 

Zerynthia r. rumina Outgroup AB095660 AF170870 1,4 











1: Caterino et al. (2001); 2: Caterino & Sperling (1999); 3: Makita et al. (1999); 4: Omoto et al. (2004); 5: 
Vila & Bjorklund (2004); 6: Yagi et al. (1999). 


(Saitou and Nei, 1987) was performed using a Kimura 2-parameter model (Kimura, 1980). 
Confidence in each node was assessed by bootstrapping (1000 replicates). 


Results and discussion 


Genetic status of L. c. huashanensis 


A phylogenetic tree based on the NDS gene sequences (Fig. 2) is consistent with that previ- 
ously reported (Makita et al., 2000). The sequences of L. c. huashanensis possesse four nu- 
cleotide differences from that of L. c. chinensis and three nucleotide differences from that of 
L. longicaudata. 


A phylogenetic tree based on the COI gene sequences (Fig. 3) showed that Luehdorfia 
japonica Leech and L. puziloi Erschoff are included in one group and that L. longicaudata 
and L. chinensis group form another group. Furthermore, the tree indicated that L. c. 
huashanensis is genetically closer to L. longicaudata than to L. c. chinensis. The sequences 
of L. c. huashanensis possess six nucleotide differences from that of L. c. chinensis and 
seven nucleotide difference from that of L. longicaudata. 


The genetic divergence among the L. longicaudata, L. c. huashanensis and L. c. chinensis is 
supported by higher (>60%) bootstrap assay values in both NDS and COI. On the other 
hand, the bootstrapping assay value between L. p. yessoensis and L. p. inexpecta are so 
small that the divergence is unclear. The NJ tree based on the combined COI and NDS se- 
quence data (846 nucleotides) indicated that L. c. huashanensis and L. longicaudata are 
shown to be monophyletic (Fig. 4), indicating smilar topology as in the COJ tree. 


It is often difficult to determine which phylogenetic methods should be adopted. Recently, 
there have been attempts to analyze the phylogenetic relationship among various genera in 
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Fig. 2. Molecular phylogenetic tree of Luehdorfia based on 389 base pair sequence analysis of the 
NDS gene of mtDNA. 


Lepidoptera by mtDNA and there are some discrepancies in the results and phylogenetic 
trees depending on the genes used (Makita et al., 2000). When the molecular phylogenetic 
relationship is discussed, it is important to compare many phylogenetic trees obtained from 
various data such as morphology, isozyme, DNA sequence, etc. 


In spite of some discrepancies between the NDS tree and the COI tree, the present study 
clearly indicated that L. c. huashanensis and L. c. chinensis are genetically differentiated 
groups. 


A food plant switch is proposed to be an important factor in studying the differentiation of 
butterflies (Shir6zu, 1955). One of the factors in butterfly differentiation is food plant 
switching by gaining new character such as new metabolic pathway of repellent in the new 
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Fig. 3. Molecular phylogenetic tree of Luehdorfia based on 457 base pair sequence analysis of the 
COI gene of mtDNA. 


food plant. Shinkawa (1999) suggested that ancestral species of Luehdorfia had fed on aio- 
phyllous Heterotropa plant and differentiated into L. puziloi due to a change in its food 
plant to deciduous A. sieboldii. 


All four species of the genus Luehdorfia can be bred using A. sieboldii (Nose, 1990b). This 
suggests that the ancestral species of genus Luehdorfia, might have been feeding on decidu- 
ous A. sieboldii. L. longicaudata might have differentiated from the L. longicaudata-L. chi- 
nensis group after a switch in its food plant to deciduous Saruma henryi during the group 
differentiating from their ancestral species. Subsequently, L. c. chinensis differentiated from 
a part of the L. chinensis group due to a change in its food plant to aiophyllous A. forbesii. 
Finally, only L. c. huashanensis, which had retained the original feeding habitat, remained 
in the Qinling Mountains. Therefore, it could be speculated that L. c. huashanensis has re- 
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Fig. 4. Molecular phylogenetic tree of Luehdorfia based on combined NDS and COI sequence data 
(846 nucleotides). 


tained the oldest characteristics of the ancestral genus, Luehdorfia, such as food habit. The 
distribution of the L. chinensis group spreads along the Chang Jiang river and its tributary 
area, and there is a large nondistributed area between the habitats of the L. c. huashanensis 
and L. c. chinensis. It might be possible to hypothesize from this distribution pattern that 
the habitat of the L. chinensis group was originated from the Qinling Mountains and spread 
to the downstream area and that the big Chang Jiang river played some role in their evolu- 
tion. A few species of the ancestral genus, Luehdorfia, have spread their habitats to the 
north while retaining their original feeding habitat and were differentiated into the L. japon- 
ica-L. puziloi group. Later, a group that changed their food plant to aiophyllous Hetero- 
tropa species was L. japonica. 
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Fig. 5. Habitat of Luehdorfia chinensis huashanensis. 


Conservation of L. c. huashanensis 


The International Union for Conservation of Nature and Natural Resources (IUCN) Red 
Data Book (Collins & Morris, 1985) annotated Luehdorfia chinensis as a data deficient 
species. Following various ecological studies (Lee, 1978; Ohno & Matsunaga, 1981; Nose, 
1990a; Hu et al. 1992; Watanabe, 1996), Yuan et al. (1998) argued that this species should 
be classified as a vulnerable species, based on their evaluations according to criteria Ala, c, 
d, and A2c of the IUCN category for endangered species (IUCN, 1994). The Chinese gov- 
ernment designated L. c. huashanensis as a class II key state protection animal in the List of 
Critical Wild Animals for Conservation based on the Wild Animal Conservation Law of the 
People’s Republic of China, enforced on March 1, 1989. In China, recent environmental 
changes in the habitat of L. chinensis have led to habitat loss and degradation and over-har- 
vesting of its host plant for medicinal use poses a threat to the survival of the species (Yuan 
et al., 1998). With regard to the conservation of the family Papilionidae, the IUCN stated 
that “it may be unwise at this stage to divert extremely restricted resources from the conser- 
vation of species to conservation of subspecies ... undue emphasis on these is not at present 
feasible, unless local expertise and enthusiasm which cannot be channeled to higher priori- 
ties are available” (New & Collins, 1991). The IUCN also emphasized the extreme impor- 
tance of accumulating biological data and the need for enthusiasm from the local re- 
searchers for conservation at the subspecies level. 


It is very important to clarify the phylogenetic specificity of subspecies huashanensis to en- 
sure its conservation in the native habitat, i. e. the Qinling Mountains, where L. c. huasha- 
nensis is locally distributed. Butterflies of this subspecies reside in the remaining patches 
of deciduous forest habitat surrounded by afforested coniferous forests (Fig. 5). 
Additionally, their food plant, A. sieboldii, is collected for use in Chinese herbal medicine, 
and there are now areas in the mountains where it is difficult to even find the butterflies 
(Yuan et al., 1998). It is known that local people are over-collecting these butterflies from 
some habitats as a source of income (Yuan et al., 1998). To prevent extinction of this rare 
subspecies, protective measures need to be implemented immediately. A few conceivable 
countermeasures are as follows: 1) preservation of habitats, 2) increasing food plants, and 
3) education of local people with respect to the conservation of nature. 
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We expect that the current molecular phylogenetic study of mtDNA will bring attention to 
the importance of L. c. huashanensis and that it can actually initiate conservation work. 
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